Two component dynamics of the superconducting order parameter revealed by 

time-resolved Raman scattering 
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We study the dynamics of the superconducting order parameter in the high-T c cuprate 
Bi 2 Sr2CaCu208-(5 by employing a novel time-resolved pump-probe Raman experiment. We find 
two different coupling mechanisms that contribute equally to the pair breaking peak. One coupling 
sets in very fast at 2 ps and relaxes slow, while the other one is delayed and sets in roughly at 5 ps 
and relaxes fast. A model that couples holes through phonons is able to reproduce one part of the 
condensate dynamics, thus, we argue that hole-spin interactions are of importance as well. 
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The nature of the interaction between holes leading to 
superconductivity is encoded in the properties of the su- 
perconducting order parameter [TJ [2] . These properties 
are reflected by the energy, the momentum dependence, 
and the time scales on which the order parameter reacts 
to an external perturbation [1 ; 2, 3,0]. In a material with 
competing interactions, there is a potential for the devel- 
opment of competing ordering phenomena [2j [3]. Un- 
doped high temperature superconductors are antiferro- 
magnetic insulators that become superconducting upon 
doping. In the superconducting state, the suppressed an- 
tiferromagnetic order might prevail on short length and 
time scales and, hence, affects thermodynamic proper- 
ties [U El [3]. Therefore, it is crucial to understand the 
transient physics that is directly connected to a specific 
phase transition [4] . Such an approach allows the study of 
competing order parameters, their individual relaxation 
channels and elucidates the potential interplay between 
them. 

Major progress in the development of pulsed laser, syn- 
chrotron, and free electron laser sources has led to in- 
novative time-resolved techniques that can directly ad- 
dress the transient physics of, e.g. a correlated mate- 
rial [5j [6j [TJ EJ El US]- However, several pump-probe 
techniques that have been applied to understand the su- 
perconducting condensate in high-temperature supercon- 
ductors, fail to probe the superconducting order param- 
eter A(k) directly [6 . In this context an intense debate 
exists on the question of whether or not one deals with 
one or more coupling processes for the pairing mecha- 
nism [9 . The two mechanisms which attracted the great- 
est attention invoke the coupling between holes through 
phonons or spin fluctuations [TJ O [3l [IT]. However, un- 
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til today there is no clear understanding to which degree 
these mechanisms jointly contribute to the superconduct- 
ing state. Different coupling mechanisms can be expected 
in different time scales in the response to an external per- 
turbation. 

In this letter, we present a unique time-resolved 
two-color inelastic light (Raman) scattering experiment, 
which allows us to probe directly the superconducting 
order parameter. We report measurements of the dy- 
namics of A(k). The samples are slightly overdoped 
Bi2Sr2CaCu208-5 (Bi-2212) high-temperature supercon- 
ductors with a T c of 82 K and are very well characterized 
Q21 [131 HH [15] . In our pump-probe Raman experiment 
we employ the UT-3 Raman spectrometer, which is a 
fully reflective achromatic spectrometer for the frequency 
range from deep ultraviolet to near infrared [16 . In order 
to obtain well defined time resolution we use two pulsed 
laser beams. The pump beam at 3.44 eV photon energy 
[full width at half maximum (FWHM) is 0.8 ps] initiates 
the experiment and drives the system out of its equilib- 
rium state. With the Raman probe beam at 1.72 eV pho- 
ton energy (FWHM = 0.9 ps) we observe the energy and 
height of the pair breaking peak A(k) as a function of de- 
lay, i.e. time difference between the pump and the probe 
beam. The large energy difference between the pump and 
the probe beam avoids any spurious signal in the Raman 
probe. We have employed Bi^-polarization ex (x 2 -y 2 ) re- 
sembling the d-wave symmetry of the superconducting 
order parameter. B^-symmetry can be studied by using 
crossed polarization between incident and scattered light 
with respect to the a and b axes in the Cu02 planes. The 
pump beam heats the sample roughly 40-65 K above its 
equilibrium temperature, whereas the probe beam heats 
the sample by about 3 K. These values are estimated 
by comparing the non-equilibrium Raman spectra to the 
corresponding temperature corrected equilibrium Raman 
spectra from a continuous wave (CW) source, and numer- 
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FIG. 1: (color online) Steady state and time-resolved Raman 
spectra of superconductors and metals depending on temper- 
ature in Bi^ geometry, (a) shows the spontaneous Raman 
scattering intensity at 10 K and at 300 K. A gap opens be- 
low 250 cm -1 (blue area) and a pair breaking peak appears 
around 420 cm -1 (red area), (b) Temperature evolution of the 
Raman difference spectrum at a delay between the pump and 
the probe beam of 3 ps. Positive values indicate an increase 
of intensity in the pumped state compared to the equilib- 
rium state, negative values indicate the opposite. The inset 
shows the integrated intensity of the difference spectra be- 
tween 300 cm -1 and 600 cm -1 . 
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FIG. 2: (color online) Temporal evolution of the time-resolved 
Raman difference spectra at an equilibrium temperature of 
10 K in Bi^ geometry. In (a) three difference spectra of three 
different delay times are shown. The contour plot shown in 
(b) presents 12 Raman difference spectra for different delay 
times. The dashed line separates two energy regions of the 
pair-breaking peak that reveal different characteristic behav- 
ior. The intensity changes are color coded demonstrating the 
transfer of spectral weight from high to low energies after 1 ps, 
respectively. 



ical estimates [TTJUH]. In the superconducting state the 
temperature rise is large enough to break Cooper pairs. 
Thus, changes in the energy and the height of the pair- 
breaking peak, i.e. the superconducting order parameter, 
can be probed. 

In Fig. [IJa) we show Bose-corrected steady state 
spectra for a transferred energy (Raman shift) between 
100 cm -1 and 600 cm -1 in the normal and superconduct- 
ing state of Bi-2212. We employ a probe energy of 1.72 eV 
in order to minimize resonance effects that are known to 
occur for higher incident photon energies [I2j [19] . In the 
normal state one can observe a flat background (black 
curve) that results from the scattering of charges within 
a marginal Fermi liquid p~9j [20j [2TJ [22] . On the other 
hand, in the superconducting state a gap opens (blue 
shaded area) and a pair breaking peak forms at roughly 
twice the maximum value of A(k) (red shaded area) at « 
420 cm- 1 ^ 52 meV [19 . The height of the peak is pro- 
portional to the number of Cooper pairs that are broken 
around (±7r,0) and (0,±7r) in the Brillouin zone. 

If we pump the superconducting state, we raise the 
sample temperature. For different equilibrium tempera- 
tures and at a fixed delay of 3 ps Raman difference spec- 
tra normalized to the integrated scattering intensity Ii nt 
are shown in Fig. [ljb) . In the pumped state, we clearly 
find a difference spectrum with positive values at low 
energies and negative values at higher energies. This 
observation reflects directly a loss of Cooper pairs and 
correspondingly a decreased intensity in the pair break- 
ing peak. Due to charge conservation, this decrease of 
intensity is accompanied by an increase of quasi-particle 
spectral weight within the gap. Finally, the Raman dif- 
ference spectrum vanishes at T c . The integral over the 
magnitude of the Raman difference spectrum is shown 



in the inset of Fig. [TJb) clearly demonstrating that we 
observe effects that are related exclusively to the super- 
conducting state and are not related to the pseudogap 

001. 

In order to analyze the dynamics of Cooper pairs in 
the pumped spectrum, we show of example in Fig. [2|a) 
three Raman difference spectra at 10 K as a function of 
delay time between pump and probe pulse. Already at 
1.65 ps after the pump beam, the Raman probe detects a 
pronounced loss of spectral weight at high energies which 
corresponds to the high-energy tail of the pair-breaking 
peak shown in Fig.JTJa), and only a marginal gain of spec- 
tral weight within the gap. At 6.6 ps the loss of spectral 
weight is shifted to energies around the maximum inten- 
sity in the pair-breaking peak and a clear gain of spectral 
weight within the gap [blue shaded area in Fig. [IJa)] is 
observed. Furthermore, the integrated change at 6.6 ps is 
significantly larger as compared to the change at 1.65 ps. 
This clearly demonstrates that the superconducting con- 
densate reacts on two different time scales to the pump 
beam. Interestingly, these time scales are roughly con- 
nected to two sides of the pair-breaking peak, i.e. above 
and below 420 cm -1 . After 16.5 ps the difference spec- 
trum is overall strongly diminished indicating a complete 
relaxation back into the equilibrium state. 

In Fig. [2^b) we show the temporal evolution in more 
detail by employing a density plot consisting of 12 Ra- 
man difference spectra (green color corresponds to no 
intensity change, i.e. the equilibrium state). The Raman 
shift is set on the y-axis, whereas the time delay between 
pump and probe beams is set on the x-axis. From this 
data set we can derive the two different time scales in 
more detail: First, a fast response at around 2 ps that 
starts with a suppression of a pair breaking peak above 
420 cm -1 . This time scale relates directly to the results 
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of novel time-resolved ARPES experiments by Perfetti et 
al showing unambiguously that the hot electrons and 
hot phonons created by a pump pulse thermalize within 
50 fs and 2 ps, respectively [23]. From this we can con- 
clude that our observed dynamics after 2 ps is driven 
by a homogeneously heated sample. After the initial 
loss of spectral weight in the pair-breaking peak around 
2 ps, an increase of low-energy spectral weight is observed 
roughly 1 ps later, reflecting the transformation of holes 
from the superconducting condensate into quasiparticles 
due to charge conservation. Then, a second, delayed re- 
sponse after 5 ps is observed in the pair breaking peak, 
this time below 420 cm -1 . This spectral weight suppres- 
sion in the pair-breaking peak also yields an additional 
gain in spectral weight within the gap roughly after 1 ps. 
Thus, summarizing our spectra, we can clearly identify 
two contributions to the suppression of the pair-breaking 
peak that have their typical energy and time scales. We 
observe a fast onset of suppression of the pair-breaking 
peak above 420 cm -1 and a delayed suppression below 
420 cm -1 . Both regions of the pair-breaking peak are 
indicated by the differently red shaded areas in Fig.JTJa). 
As expected from charge conservation, both responses 
yield their respective gain of spectral weight within the 
gap indicating a clear redistribution of the superconduct- 
ing condensate [blue shaded area in Fig. [TJa)]. 

Having identified two different time scales of the super- 
conducting condensate, what is their corresponding time 
(decay) constant? To study this, we have integrated the 
spectral weight along the energy axis above and below 
420 cm -1 as indicated by the dashed line in Fig. [2^b). 
The temporal evolution of these spectral-weight changes 
are displayed in Fig. [3^ a) and (b), respectively. As we will 
argue below, it is remarkable that the behavior of the fast 
high energy response is quite consistent with a relaxation 
process through in-plane phonons. Our calculation will 
show that this characteristic time scale of the relaxation 
process of about 7.4 ps is consistent with previous esti- 
mates from other experiments [6j[9|. However, the sec- 
ond, delayed response below 420 cm -1 is unexpected and 
has also an untypical behavior, since rise and decay times 
are roughly equal and of only 1.4 ps [see Fig.[3jb)]. From 
neutron scattering experiments we know that for energies 
below 420 cm -1 the coherence peak in the spin suscepti- 
bility yields to a strong coupling between holes even close 
to T c [24 . Thus, it is reasonable to assume that remnant 
spin-mediated coupling between holes yields an enhanced 
stiffness of the paired holes against an external pertur- 
bation resulting in this delayed reduction of the Cooper 
pairs after 5 ps. Such a scenario could be supported by a 
picture in which charge rich domains are surrounded by 
fluctuating charge poor domains. Indeed, static domains 
have been observed by static probes such as Scanning 
Tunneling Microscopy (STM) experiments [25] . 

Finally, we show in Fig. a comparison between a 
model and our Raman difference spectra in the pumped 
state. For our calculations, we have employed the method 
of density-matrix theory [26 , 27 which has been recently 




Time (ps) Raman shift (cm -1 ) 



FIG. 3: (color online) Integrated intensity for the two differ- 
ent energy regions and comparison between theoretical and 
experimental response, (a) Energy region from 420 cm -1 to 
580 cm -1 , (b) Energy region from 300 cm -1 to 410 cm -1 . 
The green characteristic decay time of (a) r = 7.4 ps and (b) 
r = 1.4 ps are determined by an exponential and Gaussian 
fit, respectively, (c) Calculated time-resolved Raman differ- 
ence spectra for two different delay times. The inset shows 
the experimental difference spectra. 



generalized for time-resolved spectroscopy on supercon- 
ductors. We have used the Hamiltonian of Ref. [18] 
that invokes the coupling to the most important phonon 
modes of the copper-oxygen planes, i.e. the breathing 
and buckling modes. The corresponding electron-phonon 
matrix elements are treated within LDA. Coupling of 
holes to spin fluctuations is not considered. Then, in 
the superconducting state, the non-equilibrium Raman 
intensity can be calculated from the imaginary part of 
the response function 

X (q = 0,„) ~ ^7k(«k(4«k) + ^(l-(/3^k)) 

k 

+M((44) + (^ k ))) . (i) 

q and uo denote the transferred momentum and energy, 
respectively. The summation runs over all (electronic) 
wavevector k in the first Brillouin zone. 7k oc (k 2 — k 2 ) 
is the non-resonant B\ g Raman vertex. Due to the em- 
ployed Bogoliubov transformation the Raman suscepti- 
bility x(q = 0, uo) includes the mixing amplitudes 
and the creation and annihilation operators «k, o^, Ac 
and of the Bogoliubov quasiparticles (linear combina- 
tion of holes and electrons), respectively [28]. The use 
of the low photon energies in the probe beam of about 
1.7 eV makes this non-resonant vertex more applicable 
as compared to visible and UV-photon energies where 
strong resonance effects need to be considered in Bi-2212 
[T2] , We chose a one band tight-binding fit to the mea- 
sured band structure p~9j [29] and a d-wave order param- 
eter Ak = Aq (cosk x - cosk y )/2, yielding a quasiparti- 
cle dispersion = y (e£ — fi) 2 + A^.. Then, we calcu- 
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late the time- dependent expectation values (a^ak)^)? 

(4Ac)W and <Aca:k>W = (olPl)*(t) with the help 
of coupled Boltzmann-equations for the non-equilibrium 
situation in the pumped state. An additional damping 
S = 5 meV is used to account for other scattering pro- 
cesses that we do not take into account. After numerical 
solution of the equations of motion, the results can be 
inserted into Eq. 0, and a pump-probe difference Ra- 
man spectrum as a function of delay time can be readily 
calculated. 

Figure [3jc) shows the calculated difference Raman 
spectra at a fixed delay time of 3.3 ps and 15 ps. The 
inset shows the corresponding measured data that have 
been corrected for the ratio of the pumped to the probed 
volume [30]. We find a fair agreement considering that 
we have used no adjustable fitting parameters in the nor- 
malized plot of the Raman response. However, the relax- 
ation in the measured Raman spectra is faster as com- 
pared to the model calculations indicating a superposi- 
tion of different phonon modes that contribute to the re- 
laxation process. The electron-phonon coupling strength 
for the in-plane oxygen breathing and out-of-plane buck- 
ling modes [31] yields typical time constants of 4 ps and 
20 ps, respectively (not shown). From this we derive 
a roughly 80% contribution of the breathing mode in 
the phonon mediated relaxation process p~8j [27]. Fur- 
thermore, it is obvious that the second response below 
52 meV with a delayed onset after 5 ps and very fast de- 
cay time of less then 1.4 ps cannot be described within 
our model. As mentioned above, this strongly suggest 
that the hole-spin interaction and the inhomogeneous dis- 
tribution of holes within the copper-oxygen plane needs 



to be included in any theory that aims to fully understand 
the time-resolved Raman response as shown in Fig.[2jb). 

In conclusion, we present a unique two-color Raman 
experiment revealing the ultrafast dynamics of the su- 
perconducting order parameter in Bi-2212 by employing 
a novel time-resolved pump-probe Raman experiment. 
Our results clearly demonstrate that the pair-breaking 
peak in the Raman responses reacts on two different time 
scales. These time scales are equivalent to two different 
coupling mechanisms. Both couplings show the redistri- 
bution of spectral weight from the pair-breaking peak to 
a quasiparticle response within the gap. The first, fast 
response sets in at 2 ps and relaxes within 7.4 ps, the 
second response sets in at 5 ps and relaxes within 1.4 ps, 
respectively [see Fig.[2|b)]. Based on density- matrix cal- 
culations, we are able to model the first response by the 
coupling of holes to phonons and find the in-plane breath- 
ing mode to be dominant. However, for an understand- 
ing of the second, delayed Raman response one would 
need to account explicitly for hole-spin interactions and 
for the inhomogeneous nature of the charge distribution. 
Our new technique provides direct information on the 
condensate dynamics as well as on the coupling mecha- 
nisms responsible for the formation of the superconduct- 
ing state. 
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